PROCEEDINGS 

OF SCIENCE 



o 

O 

(N 

> 

O 

O 



> 

00 

o 



O 



X 



A, Ai and in three-flavor (lattice) QCD 



Elizabeth D. Freeland 

School of the Art Institute of Chicago, Chicago, Illinois 60603, USA 
E-mail: |eliz@fnal ■ gov 



Andreas S. Kronfeld, James N. Simone, and Ruth S. Van de Water 

Fermi National Accelerator Laboratory, Batavia, Illinois, 60510, USA 

Fermilab Lattice Collaboration with the MILC Collaboration 

The heavy-quark expansion for inclusive semi-leptonic B decays introduces A and X\ , which are 
matrix elements in heavy-quark effective field theory. We review how they can be obtained from 
an analysis of the heavy quark mass dependence of heavy-light meson masses in lattice QCD. We 
present preliminary results for the bottom quark mass, m^, using A and X\ for the B s meson from 
a 2+1 sea- flavor unquenched calculation. 



XXIVth International Symposium on Lattice Field Theory 
July 23-28, 2006 
Tucson, Arizona, USA 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



A, Ai and my in three-flavor (lattice) QCD 



Elizabeth D. Freeland 



1. Background and Overview 

The study of heavy-quark physics has progressed greatly since the discovery of heavy-quark 
spin-flavor symmetry and the development of heavy-quark effective theory (HQET) to systemat- 
ically deal with symmetry-breaking effects. HQET is based on a scale separation between the 
physics of heavy and light quarks. It yields expressions for observables as expansions in in- 
verse powers of the heavy-quark mass or, alternatively, inverse powers of the heavy-meson mass. 
These expansions share a common set of HQET matrix elements which must be evaluated non- 
perturbatively 

One would like to determine these matrix elements on general principles and, more impor- 
tantly, because they are needed to ascertain the Cabibbo-Kobayashi-Maskawa (CKM) matrix ele- 
ments via inclusive decay measurements [Q, |]]. For example, the heavy-quark expansion for the 
rate of the decay B —> X c Iv is given by [|5j] 
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where \V c b\ is the CKM matrix element of interest; Mb is the B-meson mass; and, A, k\, and 
%l are scheme-dependent hadronic matrix elements defined in HQET. Currently, the HQET ma- 
trix elements are determined by fitting measurements of various moments of heavy-meson decay 
distributions to corresponding HQET expressions [0, H,^J,|7|]. 

In addition to using experimental measurements, one would like to calculate the HQET matrix 
elements from first principles. In fact, because the same non-perturbative quantities appear in the 
HQET expression for the meson mass, there exists a rather direct method for calculating them 
using lattice QCD. This method was first proposed in Ref. D9Q, where the corresponding quenched 
calculation was reported. 

The HQET expression for the mass of a heavy-light meson is [Q, [|, ||] 

M = m + A-^-d/-^ + 0(l/m 2 ), (1.2) 
2m 2m 

where / is the total meson angular momentum, and do = 3 and d\ = — 1 for the pseudoscalar and 
vector mesons respectively. The mass of the heavy-light meson is M and that of the heavy quark is 
m. Working with the spin-averaged meson mass, M, the equation simplifies to 

M-m = A--^- + 0(\/m 2 ). (1.3) 
2m 

Since a lattice calculation allows one to work at any heavy-quark mass desired, we can generate 
data for a variety of heavy-light meson masses and then fit to Eq. ( \l.3\ > to determine A and X\ . 

The key is that the formalism of HQET applies to any underlying theory with the heavy-quark 
spin-flavor symmetry, such as a lattice gauge theory with heavy Wilson quarks [JlO , 11]. The lattice 



introduces another short distance, which can be treated perturbatively via the Wilson coefficients. 
The resulting expression for the spin-averaged meson mass on the lattice is Jl0| ] 

Mj -mj =A(a)-^^- + 0(l/m 2 ). (1.4) 
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Here, the heavy-quark a dependence is absorbed into the rest mass mi and the kinetic mass wi2, 
defined via the heavy-quark propagator [|l2|]. Discretization effects of the gluons and light quarks 
lead to a separate a dependence of the HQET matrix elements A and X\ . My is the spin-averaged 
rest mass on the lattice. 

To arrive at A and X\ for B^" and we would need an expression for the (spin-averaged) 
meson mass from chiral perturbation theory. Continuum HQET expressions exist in the litera- 



ture Q13 ]. The full expression including effects from staggered quarks and HQET is being derived 
now. Because we can do simulations with light valence quarks near or at the strange-quark mass, 
and because we find that the effect of sea quarks is mild, we obtain preliminary results for the 
bottom quark mass, rrib, using A and X\ from the B s meson. 

In the following, we first discuss the lattice calculation of meson and quark masses. We then 
discuss how M\ — m\ (and hence A and X\) depends on sea and valence quark masses and the 
lattice spacing. Finally, we calculate the binding energy for a spin-averaged B s meson and use that 
to arrive at a preliminary value for the bottom-quark mass. 

2. Meson Masses, Quark Masses: Determining M, mi and m,2 



We use the MILC unquenched gauge configurations [ ]14[ ] with 2+1 flavors of sea quarks and a 
Symanzik-improved gluon action. We use three lattice spacings: a = 0.18,0.15, and 0.12 fm. Both 
sea and light valence quarks use the "asqtad" staggered-fermion action JIB]]. Light valence quarks 
have masses ranging from m q = \ .\m s to 0.lm s , where m s is the (physical) strange quark mass. 
Masses of the two light sea quarks range from approximately 0.05 m s to 0.1 m s . For heavy quarks, 
we use the Fermilab action In anticipation of the full calculation of A and X\, we use seven 
or more heavy-quark masses at each lattice spacing. They range in mass from heavier-than-bottom 
to lighter-than-charm. 

Pseudoscalar and vector meson masses are obtained from two-point correlation function fits 
done using multi-state, constrained curve fitting [17]. Both % 2 and fit stability are used to determine 



the goodness of fit. The results are spin-averaged to obtain M\. 



Equation ( L4 ) on its own does not specify the renormalization scheme for the masses and, 



hence, A and X\ . Although it is straightforward to obtain the pole rest and kinetic masses at the one- 



loop level [|12j |18| |, the perturbative expansion of the pole mass is marred by infrared effects [|19|]. It 
is, therefore, better to introduce a short-distance mass. Because the MS mass does not run correctly 
for renormalization scales below the heavy-quark mass scale, it is not appropriate. Several other 
short-distance mass definitions are available in the literature. Here, we use the potential-subtracted 
mass, mps [pCQ, which is based on the static quark potential and introduces a separation scale, }Xf, 
where A QCD < \x f < 2GeV. 

For a s , we use the V-scheme; scale setting, q* , is done via the Brodsky-Lepage-Mackenzie 
(BLM) prescription The value of a s (q*) is obtained from the average value of the plaquette 
and the four-loop j3 -function as described in [E20. 



3. Sea Quark, Valence Quark and Lattice Spacing Dependencies 

The value of the meson binding energy, M\ —ni\ ps, depends upon the sea quark masses, the 
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Figure 1: (a) M\ — m\ps vs 1 /2ni2,ps for three values of light sea quark masses; the valence quark mass 
m q = m s and a = 0.125 fm. (b) M\ — m^ps vs l/2m.2,ps for several valence quark masses, m q . Here, 
a = 0.125 fm and m u c i = 0.18m s . 

Bottom and charm quarks have values of 1 /2rri2ps = 0. 13 and 0.58GeV _1 , respectively. 



light valence quark masses (m q ) and the lattice spacing. We will discuss each of these dependencies 
in turn and estimate numerically how they affect our results for the B s system. 

We begin with the effects of the 2+ 1 sea quarks. Figure |](a) is a plot of M\ —m\ t ps vs l/2»i2,PS 
for three values of sea-quark mass ratios m u ^/m s . For each graph, A is the intercept and X\ is the 
slope, while the curvature is related to a combination of HQET matrix elements at 0(1 /m 2 ) [^]. 
Figure [j](a) allows one to view the dependence of these quantities on the sea quarks. One can 
see that varying the light sea-quark mass has only a small effect on M\ — m\ps (or A and X\). In 
evaluating M\ — m\fs for B s we used, at each value of a, the ensemble with the lightest available 
m U}C j sea quarks, and used the variation from different ensembles in our estimate of the systematic 
error. 

Figure |I|(b) is a plot of M\ — m\ ps vs Xjlm^^ps for several light valence quark masses. As 
expected, the meson mass depends strongly on the value of the light valence quark mass. For the 
numerical results reported below, we will consider only mesons with a strange valence quark. To 
arrive at M\ — m\ ps for B s , we use the m q = 1.06ra v result and allow an uncertainty based on the 
value at m q = 0.77 'm s . 

Figure || is a plot of Mi —m^ps vs \/2ni2,ps for three lattice spacings: 0.125,0.15 and 0.1 8 fm. 
Results for the coarsest lattice spacing fall between the two finer spacings; the error bars shown, 
however, are statistical only. Discretization errors appear from two sources. The first is from the 
light quarks and gluons. These errors first appear at 0(a s a 2 ), 0(a 4 ). The second is from the trun- 
cation of the perturbative series for the heavy-quark masses mi and m%. Because the origins of the 
discretization effects from each sector are isolated, we can analyze their contributions separately. 
Nevertheless, over the range of masses and lattice spacings we are working with, the two uncer- 
tainties are of comparable size, so the behavior of the total discretization error is non-trivial. It is 
expected that the inclusion of two-loop effects would clarify the results shown here. To estimate 
the error due to truncation, we take the 2-loop contribution to be a factor of a smaller than that of 
the 1-loop contribution. Results from extrapolations of M\ — m\ ps in a 2 and a 4 were used to assign 
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Figure 2: M\ — m\ps vs 1 /2ni2,ps for three lattice spacings: 0.125,0.15 and 0.18 fm. For all cases m q = m s . 
Error bars are statistical only. Bottom and charm quarks have values of 1 /2ni2,ps = 0.13 and 0.58 GeV -1 
respectively. 



a systematic error for the discretization of light quarks and gluons. Our final error budget includes 
an uncertainty in the determination of a itself. 



4. Result for the &-quark mass 

Focusing on we calculate the binding energy of a (spin-averaged) B s meson, by averaging 
the results at the a = 0.125 fm and 0.15 fm spacings. Our result is 0.99(18) GeV in the potential- 
subtracted scheme with a factorization scale of jif = 2.0 GeV. Using this value of the binding 
energy, we can make a preliminary estimate of the value of the bottom-quark mass. 

m b = M exp - (Mi - m hPS ) (4. 1) 

where M exp is the experimentally measured, and spin-averaged, value of the B s mass; we use the 



value of 5.402(2) GeV from values from the Particle Data Group []23|]. This yields a preliminary 



value mbps = 4.41(18) GeV. We do not quote a result in the MS scheme at this time. For com- 



parison, a QCD sum rule calculation [24] obtains, m^ps = 4.52(6) GeV at \if = 2.0 GeV and 
^ms= 4 -!9(6) GeV. 

Table [T] provides a preliminary budget of the uncertainties in this calculation. The two largest 
are the uncertainty due to truncation of the QCD perturbation theory for the quark masses and the 
uncertainty due to the continuum extrapolation. All uncertainties are added in quadrature to arrive 
at the total. 



5. Summary and Outlook 

We report a preliminary calculation of the bottom quark mass using a (lattice) HQET calcu- 
lation of the spin-averaged B s binding energy. We use lattice QCD with 2+1 flavors of staggered 
sea quarks. Heavy-light mesons are constructed from a staggered valence and Fermilab heavy 
quark. We find m^ps = 4.41(18) GeV in the potential-subtracted scheme with a factorization scale 
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Source GeV 

statistical 0.005 

inputs (a, JQ,, JQrit, u$) 0.041 

sea-quark mass dependence 0.04 

strange-quark mass tuning 0.025 
perturbation theory (heavy quark discretization) 0.10 

light quark and gluon discretization 0. 14 

total 0.18 



Table 1: Uncertainties in the quantity M\ —ni\_ PS for a spin-averaged B s meson. 

of jU/ = 2.0 GeV. The dominant uncertainties in this calculation can be reduced by the inclusion 
of 2-loop effects in the perturbative expansions for m\ and mi, and with improved understanding 
of light quark and gluon discretization effects. Future work will include the calculation of HQET 
matrix elements A and X\ for B ± and Bj, which can be used in the determination of CKM matrix 
elements from their inclusive, semileptonic decays. 

Acknowledgments 

E.D.F. would like to thank the Fermilab Theoretical Physics Department for their hospitality 
while this work was being done, and Don Holmgren and Amitoj Singh for their help with comput- 
ing issues. E.D.F. is supported in part by the American Association of University Women. Fermilab 
is operated by Universities Research Association Inc., under contract with the DOE. 

References 

[1] E. Eichten, Nucl. Phys. Proc. Suppl. 4, 170 (1988); E. Eichten and B. Hill, Phys. Lett. B 243, 427 
(1990). For a pedagogical introduction see: A. F. Falk, hep-ph/0007339. 

[2] A. F. Falk and M. Neubert, Phys. Rev. D 47, 2965 (1993) [hep-ph/9209268]. 

[3] J. Chay, H. Georgi and B. Grinstein, Phys. Lett. B 247, 399 (1990); 1. 1. Y. Bigi, N. G. Uraltsev and 

A. I. Vainshtein, Phys. Lett. B 293, 430 (1992) [Erratum-ibid. B 297, 477 (1993)] [hep-ph/9207214]; 
1. 1. Y. Bigi, M. A. Shifman, N. G. Uraltsev and A. I. Vainshtein, Phys. Rev. Lett. 71, 496 (1993) 
[hep-ph/9304225]; B. Blok, L. Koyrakh, M. A. Shifman and A. I. Vainshtein, Phys. Rev. D 49, 3356 
(1994) [Erratum-ibid. D 50, 3572 (1994)] [hep-ph/9307247]; T. Mannel, Nucl. Phys. B 413, 396 
(1994); [hep-ph/9308262]. 

[4] M. Gremm and A. Kapustin, Phys. Rev. D 55, 6924 (1997) [hep-ph/9603448] . 

[5] D. Cronin-Hennessy et al. [CLEO Collaboration], Phys. Rev. Lett. 87, 251808 (2001) 
[hep-ex/0108033]. 

[6] S. Chen et al. [CLEO Collaboration], Phys. Rev. Lett. 87, 251807 (2001) [hep-ex/0 108032]; 

B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 93, 01 1803 (2004) [hep-ex/0404017]. 

[7] M. Gremm, A. Kapustin, Z. Ligeti and M. B. Wise, Phys. Rev. Lett. 77, 20 (1996) [hep-ph/9603314]; 
M. Gremm and I. Stewart, Phys. Rev. D 55, 1226 (1997) [hep-ph/9609341]; A. F. Falk and 
M. E. Luke, Phys. Rev. D 57, 424 (1998) [hep-ph/9708327]. 



6 



A, X\ and nib in three-flavor (lattice) QCD 



Elizabeth D. Freeland 



[8] T. Mannel, Phys. Rev. D 50, 428 (1994) [hep-ph/9403249]; 1. 1. Y. Bigi, M. A. Shifman, 
N. G. Uraltsev and A. I. Vainshtein, Phys. Rev. D 52, 196 (1995) [hep-ph/9405410]. 

[9] A. S. Kronfeld and J. N. Simone, Phys. Lett. B 490, 228 (2000) [Erratum-ibid. B 495, 441 (2000)] 
[hep-ph/0006345]. 

[10] A. S. Kronfeld, Phys. Rev. D 62, 014505 (2000) [hep-lat/0002008]. 

[11] J. Harada, S. Hashimoto, K. I. Ishikawa, A. S. Kronfeld, T. Onogi and N. Yamada, Phys. Rev. D 65, 
094513 (2002) [Erratum-ibid. D 71, 019903 (2005)] [hep-lat/01 12044]; J. Harada, S. Hashimoto, 
A. S. Kronfeld and T. Onogi, Phys. Rev. D 65, 094514 (2002) [hep-lat/01 12045]. 

[12] B. P. G. Mertens, A. S. Kronfeld and A. X. El-Khadra, Phys. Rev. D 58, 034505 (1998) 
[hep-lat/97 12024]. 

[13] E. Jenkins, Nucl. Phys. B 412, 181 (1994) [hep-ph/92 12295]. 

[14] C. W. Bernard et al, Phys. Rev. D 64, 054506 (2001) [hep-lat/01 04002]; C. Aubin et al, Phys. Rev. D 
70, 094505 (2004) [hep-lat/0402030]. 

[15] T. Blum et al, Phys. Rev. D 55, 1 133 (1997) [hep-lat/9609036]. K. Orginos and D. Toussaint [MILC 
collaboration], Phys. Rev. D 59, 014501 (1999) [hep-lat/9805009]. J. F. Lagae and D. K. Sinclair, 
Phys. Rev. D 59, 01451 1 (1999) [hep-lat/9806014]. G. P. Lepage, Phys. Rev. D 59, 074502 (1999) 
[hep-lat/9809157]. K. Orginos, D. Toussaint and R. L. Sugar [MILC Collaboration], Phys. Rev. D 60, 
054503 (1999) [hep-lat/9903032]. C. W. Bernard et al. [MILC Collaboration], Phys. Rev. D 61, 
1 1 1502 (2000) [hep-lat/9912018]. 

[16] A. X. El-Khadra, A. S. Kronfeld and P. B. Mackenzie, Phys. Rev. D 55, 3933 (1997) 
[hep-lat/9604004]. 

[17] G. P. Lepage, B. Clark, C. T. H. Davies, K. Hornbostel, P. B. Mackenzie, C. Morningstar and 
H. Trottier, Nucl. Phys. Proc. Suppl. 106, 12 (2002) [hep-lat/01 10175]. 

[18] M. Nobes and H. Trottier, PoS LAT2005, 209 (2006) [hep-lat/0509128] and private communication. 

[19] M. Beneke and V. M. Braun, Nucl. Phys. B 426, 301 (1994) [hep-ph/9402364]. 1. 1. Y. Bigi, 

M. A. Shifman, N. G. Uraltsev and A. I. Vainshtein, Phys. Rev. D 50, 2234 (1994) [hep-ph/9402360]. 

[20] M. Beneke, Phys. Lett. B 434, 115 (1998) [hep-ph/9804241]. 

[21] S. J. Brodsky, G. P. Lepage and P. B. Mackenzie, Phys. Rev. D 28, 228 (1983); G. P. Lepage and 
P. B. Mackenzie, Phys. Rev. D 48, 2250 (1993) [hep-lat/9209022] . 

[22] Q. Mason et al. [HPQCD Collaboration], Phys. Rev. Lett. 95, 052002 (2005) [hep-lat/0503005]; 
G. P. Lepage, private communication. 

[23] W. M. Yao et al. [Particle Data Group], J. Phys. G 33, 1 (2006). 

[24] A. Pineda and A. Signer, Phys. Rev. D 73, 1 1 1501 (2006) [hep-ph/0601 185]. 



7 



